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THREE-DIMENSIONAL  STRESSES  IN  A  HALF  SPACE 
CAUSED  BY  PENNY-SHAPED  INCLUSIONS 


INTRODUCTION 


Elastic  fields  caused  by  inclusions  in  infinite  media  have  been  extensivley  investigated  by  several 
authors  [1-5]  after  Eshelby’s  work  [6-8].  Other  research  efforts  have  addressed  the  half-space  prob¬ 
lem  with  an  inclusion  located  near  the  free  surface  [9-12],  In  these  studies,  the  following  methods 
were  used:  Galerkin  vector  [9],  Papkovich-Neuber  displacement  potential  [10],  image  stress  caused 
by  two  cuboidal  inclusions  with  uniform  eigenstrains  [11],  and  Green’s  function  in  the  half  space 
[12].  Mura  has  recently  reviewed  these  research  efforts  [13]. 

When  the  elastic  moduli  of  an  ellipsoidal  subdomain  of  a  material  differs  from  those  of  the 
remainder  (matrix),  the  subdomain  is  called  an  ellipsoidal  inhomogeneity.  Cracks,  voids,  and  precipi¬ 
tates  are  examples  of  these  inhomogeneities.  A  material  containing  inhomogeneities  is  assumed  to  be 
free  from  any  stress  field  unless  an  external  stress  field  a-j  is  applied.  On  the  other  hand,  a  material 
containing  inclusions  is  subjected  to  an  internal  stress  caused  by  the  eigenstrain  ej}  even  if  it  is  free 
from  any  external  loads.  The  definition  of  eigenstrains  has  been  given  by  Mura  [13]  and  is  the  same 
as  the  stress-free-transformation  strain  described  by  Eshelby  [6]. 

The  solutions  for  ellipsoidal  inhomogeneities  can  be  reduced  to  the  penny-shaped  or  elliptical 
crack  case  by  setting  the  elastic  constants  X  and  n  for  the  inhomogeneities  equal  to  zero.  The  solution 
of  the  three-dimensional  problems  for  these  cracks  has  received  considerable  attention  [14-19].  The 
stress  field  of  a  penny-shaped  crack  in  the  half  space  can  be  solved  by  obtaining  the  relevant  system 
of  integral  equations  for  the  problem  formulated  by  Erdogan  and  Gupta  [20]  for  the  stress  analysis  of 
multilayered  composites  with  a  flaw. 


In  the  present  study,  Eshelby’s  method  for  ellipsoidal  inclusions  [6-8]  and  Hankel’s  transforma¬ 
tion  method,  used  to  obtain  the  elastic  solutions  of  a  circular  dislocation  loop  in  an  unbounded  media 
[21]  and  in  the  half  space  [22],  are  used  for  the  analysis  of  the  elastic  solution  of  axisymmetric  inclu¬ 
sions  and  axisymmetric-ellipsoidal  inhomogeneities  in  the  half  space.  The  method  provides  a  novel 
way  for  obtaining  the  image  stresses  of  an  ellipsoidal  inclusion  in  the  half  space.  It  is  used  to  find  a 
more  general  solution  of  an  ellipsoidal  inclusion  with  anisotropic  eigenstrain.  Existing  solutions  are 
shown  to  be  special  cases  of  the  present  result.  This  method  can  also  be  used  to  obtain  the  stress 
field  of  a  penny-shaped  crack  in  the  half  space. 

BASIC  APPROACH 


In  this  report,  wc  consider  an  axisymmetric  ellipsoidal  inclusion  in  a  half  space  (Fig.  1). 
general,  the  inclusion  fi]  is  given  by 
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Symmetry  with  respect  to  the  Jt3-axis  is  then  defined  by  a,  =  a2,  and  the  anisotropic  eigenstrain 
of  the  inclusion 


ejj  =  bjjie  +  b  6, 3)  ij  =  1,  2,  3, 


(2) 


where  6,7  is  Kronecker  delta.  (Note  that  the  usual  summation  convention  does  not  apply  to  any  of  the 
expressions  in  this  report.)  Equation  (2)  states  that  only  normal  eigenstrains  appear,  and 
e] |  =  ej2  =  e  and  e33  =  e  +  b . 


For  the  inclusion  defined  by  Eq.  (1)  with  the  uniform  eigenstrain  described  by  Eq.  (2)  with 
a3  —  0,  the  stress  field  in  the  unbounded  medium  outside  ft,  is  obtained  by  using  Eshelby’s  method 
(6-8).  The  result  is  given  by 
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where  the  numerical  suffixes,  i,  j  =  1,  2,  3,  following  a  comma  denote  differentiation  with  respect 
to  the  Cartesian  coordinates  xt,x2,x3,  e.g.  4>,,j  =  d2<t>/dxjdxj,  and  <p  is  the  Newtonian  potential 
function  that  is  given  by 


2  t  V  j 
*  =  rafa  3  J  —  ds  , 


where 


*1  +  x22 
-> 

r2 

+  ,  3 

V 

a2  +  s 

a3  +  s 

A  =  (a\  +  s)(aj  +  s)1/2, 

and  Xq  is  the  largest  root  of  U  =  0  outside  of  and  Xo  =  0  inside  of  For  inclusions  with  uni¬ 
form  dilatation  eigenstrain  only  ( b  =  0),  Eq.  (3)  is  valid  for  any  a3  value.  The  detailed  expression 
of  cf>  for  both  the  oblate  spheroid  (a  j  >  a3)  and  the  prolate  spheroid  (a ,  <  a3)  are  given  by  Yu  [23]. 
Equation  (3)  can  be  transformed  into  cylindrical  coordinates  (r,  9,  z)  as  follows: 


4ir(l  —  v) 


°rr  =  -7 lfr~:  I  +  Z<t>'xz  +  ~~  ^’r  +  t  b'n 


^(1  +  v)e  , 

+  -TT, - r  -  +  <P>zz  . 

2tt(1  -  v)  r  J 


& 

4ir  (/  —  p) 


.  1  -  2P 

ZvVjzz - <P,r 
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/i(l  +  v)e 
27r(l  —  v)  r 


|K/n 

m 


—  tj.  ,  ,  u(l  +  v)e 

0tz  4tt(1  -  v)  l<i>,zz  Z*,zzz]  "  2tt(1  -  v)  *’«  ' 


„  ^  r- '  1  g(l  ±  p)g  t 

-  p)lZ*,rzz]  2tt(1  -  0,rz  ’ 


CTr#  -  -  0. 

Equations  (5)  are  obtained  with  the  aid  of  the  following  relationships: 
V24>  =  0, 

X\4>,  2  =X20,1, 


0,r  =  —  (X  X2<f>,2), 
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where  the  letter  suffixes  following  a  comma  denote  differentiation  with  respect  to  the  cylindrical  coor¬ 
dinates  r,  0,  and  z,  e.g.  <j>,n  =  d2<t>/drdz. 

For  a  circular-edge  dislocation  loop  with  the  ^-axis  as  the  axis  of  symmetry  in  an  unbounded 
medium  (Fig.  2),  the  stress  field  is  found  by  Kroupa  [21]  by  using  Hankel  transformations.  For 
z  >  0,  Kroupa’ s  solution  can  be  rewritten  as 

arr  =  ~  2(1  v)  a  '^,ZZ  +  l^,zzz  +  ~ — ““  (V')>r  +  ~(/o r'  )>rr]  - 

°ee  =  ~  ~  a  dolK  ~  f  dol)m  . 

°ZZ  =  V)  a  ~  Z^°  (7) 

an  =  2(\b-  v)  a  Iz(/0-1)»»al’ 

°re  =  °ze  =  0, 

where 

Inm  =  (rt/a)Jx{t)e-*,adt, 

C  =  -a(/rV 

=  ~arm-\r-m  +  '  InmZ\),r  (m  =  0,  1,  2. . .  ;  n  =  -1,0,  1,2...), 

and  Jm  is  the  Bessel  functin  of  the  mth  order,  a  is  the  radius  of  the  circular  dislocation  loop,  and  b' 
is  the  Burger’s  vector.  Equation  (7)  is  obtained  by  the  method  of  Hankel  transformation  as  used  for 
cylindrically  symmetric  problems  of  the  theory  of  elasticity  in  Sneddon’s  book  [24]  and  subjected  to 
the  following  boundary  conditions: 

uz(r,  0)  =  —  b '  for  0  <  r  <  a  , 

=  0  for  r  >  a  ,  (8) 

an  (r ,  0)  =  0  for  0  <  r  <  oo  . 


For  the  penny-shaped  inclusion  without  shear  and  dilatation  eigenstrains  (penny-shaped  prismatic 
inclusion),  which  is  the  axisymmetric  inclusion  when  a3  approaches  zero,  the  equivalent  eigenstrains 
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are  en  —  e 22  =  ^33  ^  0.  If  we  reduce  Eq.  (5)  for  a  penny-shaped  prismatic  inclusion,  that  is, 

a 3  —  O  and  e  =  0,  it  is  interesting  to  note  the  similarity  between  Eqs.  (5)  and  (7).  By  putting 

</>  =  */-'  (a 3  —  0) ,  (9) 

where  k  -  2irb  'a  /b ,  and  a,  =  a,  the  elastic  solutions  of  both  the  penny-shaped  prismatic  inclusion 
(Eq.  5)  and  the  circular-edge  dislocation  loop  (Eq.  7)  are  identical.  This  suggests  that  the  method 
used  to  investigate  the  elastic  solution  of  a  circular-edge  dislocation  loop  in  the  half  space  [22]  can  be 
applied  to  solve  the  elastic  field  caused  by  an  axisymmetrical  inclusion  in  the  half  space.  This 
approach  is  resonable  since  the  solution  of  the  axisymmetrical  inclusion  can  be  obtained  by  the 
integration  of  the  results  of  a  penny-shaped  prismatic  inclusion  and  the  fact  that  if  the  inclusion  has 
the  same  elastic  moduli  as  the  matrix,  the  stress  field  is  the  same  as  that  of  a  small  dislocation  loop 
when  both  the  dislocation  loop  and  the  inclusion  are  infinitesimally  small  [8].  For  example,  a  small 
inclusion  of  volume  V  and  an  eigenstrain  e\3  in  the  x3  direction  has  the  same  stress  field  as  that  of  a 
prismatic  interstitial  dislocation  loop  of  area  A  and  Burgers  vector  bx  provided  that  Ve33  =  Abx. 


Consider  the  half  space  -c3  =  z  >0  (Fig.  1),  an  axisymmetric  inclusion  with  the  center  at  the 
point  (0,  0,  c)  in  such  a  way  that  its  axis  of  symmetry  (z-axis)  is  perpendicular  to  the  plane  of  the 
free  surface  z  =  0.  In  order  that  the  plane  z  =  0  be  a  free  surface,  no  force  must  act  on  it,  thus  the 
stress  components  at  z  =  0  must  satisfy  the  boundary  conditions 


(®rz)z-  o  =  0, 


(°zz  )z  =o  =  0, 


and  the  equilibrium  condition 


£  °ij,j  =  °- 
3=1 


Y«\ 

y** 
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Similar  to  the  work  of  Bastecka  [22],  the  stress  atj  outside  the  axisymmetric  ellipsoidal  inclusion 
centered  at  the  point  (0,  0,  c)  but  in  the  half  space  z  >  0  is 


°U  =  °'j  +  +  °lj  . 


which  will  satisfy  the  required  boundary  conditions  (Eq.  10)  and  the  equilibrium  condition  (Eq.  11). 
This  converges  to  zero  for  and  x2  approaching  ±  oo  and  x3  approaching  oo.  in  Eq.  (12),  the  term 


is  the  stress  caused  by  the  axisymmetric  inclusion  (and  outside  of  it)  centered  at  the  point 


(0,  0,  c):  a'j  is  the  stress  caused  by  the  image  inclusion  Q2  centered  at  the  point  (0,  0,  -c)  (Fig.  3) 
with  eigenstrain 


-{?!])'  =  -&ij(e  +  bbn). 


l  \ 

I  (0,0, -C)  ' 


\  I  / 

\  I  / 

\  I  / 


Fig.  3  —  Semi-infinite  solid  containing  an  ellipsoidal  inclusion  fi, 
and  its  image  fi2 


VlVAs 


NRL  REPORT  9134 


Equation  13  shows  that  oj:  is  an  additional  stress  that  satisfies  the  boundary  condition 

(oL );  =o  =  -(4  +  47)z=o  =  0;  (14a) 

=0  =  "(4  +  arJz=0-  (14b) 


The  solutions  for  the  stresses  a jj  and  o-J  are  obtained  by  translating  the  origin  of  coordinates  in 
Eq.  (3)  and  Eq.  (5)  to  points  (0,  0,  c),  and  (0,  0,  —(c)  respectively.  The  Newtonian  potential  func¬ 
tion  4  and  <t>n  for  the  solutions  of  a'j  and  a-j  respectively  are  given  by 


/  2  f°°  V\ 

<t>  =  ira  fa  3  ]  -  ds 

J"i  A 


U2 


4>"  =  ira  \a  3  f - ds 

A 


(15) 


where 


U  i  =  1 


U7  =  1 


xf  +  x2 

aj  +  s 

xf  +  x\ 
a  i  +  5 


C*3  ~  f )' 
+  S 

(*3  +  C)' 


A  =  (a  ,2  +  j )  (af  +  s)1/2  , 
and  where 


X.  is  the  largest  root  of  U x  =  0  for  exterior  points  of  fl| , 


Xt  =  0  for  interior  points  of  and 


X2  is  the  largest  root  of  U2  =  0 . 


SOLUTION  FOR  o', 


Substituting  Eqs.  (5)  and  (15)  into  Eq.  (14)  gives 

(°zzh  =0  =  0: 


(°'rzh=  0  ~ 


2ir(l  —  v) 


cb<j>"  -  2(1  +  if)e<t>," 


z  =0 


(16a) 

(16b) 
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where  for  z  =  0,  </>,4;  =  and  <£,,{,  =  Now,  in  tne  limit  when  a3  approaches  zero,  that 

is,  for  the  penny-shaped  inclusion  (uj  =  a2  =  a),  we  can  substitute  Eq.  (9)  into  Eq.  (16b)  to  obtain 


(<)j  =o 


I ibc 

2ir(l  —  v) 


4  L°°  t2J \(pt )J ,(/ )e  ~c,/adt 
a  u 


where  p  =  r  /a  . 


p(l  +  v)e 
7r(l  -  v) 


4  C  tJ^pt)J{(t)e-c,/adt  , 
a  0 


(16c) 


For  the  axisymmetric  problem,  by  the  appropriate  expression  of  the  elastic  displacements  as  the 
derivatives  of  certair  function  \j/(r  ,z)  in  cylindrical  coordinates,  the  equilibrium  and  Beltrami  equa¬ 
tions  are  replaced  by  a  single  equation  [24] 


vV(r,z)  =  0,  (17) 

whose  general  solution  is  carried  out  by  the  method  of  integral  transformations.  The  function  \p  is 
replaced  by  its  Hankel  transform  of  zero  th  order, 

G(f,z)  =  |Q  r\p(r  ,z)J0(£r )  dr  ,  (18) 

and  it  can  be  shown  that  G(f,z)  is  generally  given  by  the  expression 

G(f,z)  =  (A  +  Bz)e~h  +  (C  +  Dz)e*z  ,  (19) 

where  A,  B,  C  and  P  are  unknown  functions  of  f,  which  are  determined  from  the  boundary  condi¬ 
tions.  The  stress  components  are  expressed  by  means  of  the  function  G(f,z). 

In  the  present  case,  we  consider  the  solution  to  converge  to  zero  for  z  approaching  oo.  Thus 
we  set  C  =  D  =  0.  To  determine  A  and  B  from  the  first  boundary  condition  (Eq.  16a),  we  obtain 
the  following  relationship 


A 


t  L 

\  +  H  f 


(20) 


where  X  =  2fiv/(l  -  2v)  is  Lame’s  constant.  From  the  second  boundary  condition  (Eq.  16b),  as 
modified  it;  Eq.  (16c),  we  have 


,l«. 

k<\ 

r. 


<■ 


{a'r- ).  =0  =f(r) 

=  rF(ny,(rr)dr. 


(21) 


where 


F(f)  =  -2(X  +  M)r25(n, 


(22) 
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fir)  = 


tibc  k 


2ir(l  -  ,)  a3 
Ml  3-  v)e  k 


—  J  tzJ ^t)J i(t)e  cUadt 


tJ \{pt)J  \(t)e~a  la dt . 

t(1  -  v)  a 2  0 


(23) 


By  letting  t  =  a  f,  Eqs.  (21)  and  (23)  give 

k  n 


Fit)  =  - 


2ir(l  —  v) 


[cb£  +  2(1  +  v)e]J](ane 


-cl 


(24) 


By  substiti/ing  Eqs.  (24),  (22),  and  (20)  into  Eq.  (19),  the  function  G(f,z)  is  found  that  can  then  be 
substituted  in  the  expressions  for  the  stress  a/,  [24,  §51].  After  substituting  the  relationship  again. 


4>"  =  k(l0-')“  =  k  $o  rlJfot)Jx(t)e-,lz+e)'adt  , 


(25) 


these  stresses  o' ,  are  as  follows. 
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When  e  =  0  and  4>n  =  k(I0  1  )u ,  Eq.  (26)  reduces  to  the  same  results  obtained  by  Bastecka  [22]  for 
a  circular-edge  dislocation  loop  in  the  half  space.  In  Cartesian  coordinates  Eq  (26)  becomes 


a[j  = 


fibc 


2tt(  1  -  v) 


(1  -  2v)(83i  +  8y  -  l)<f>,jj3  -  0$ 


+  2v8jj<t>,333  -  x3<t>,[j33 


fi(  1  4-  v)e 
*0  -  v) 


[(1  -  2v)(b3i  +  63j 


1  ~  <t>JJ 


+  2v6ij<t>,{[  -  x3<t>,lj3  . 

It  can  be  shown  that  a-,  satisfies  the  equation  of  equilibrium,  that  is, 

b  °!jj  =  0. 

7  =  1 


(27) 


(28) 


Therefore,  for  points  outside  Qj,  the  stress  field  caused  by  the  presence  of  a  penny-shaped  inclu¬ 
sion  in  the  half  space  can  be  obtained  by  Eqs.  (3),  (12),  and  (27).  Thus, 

Oij  =  4T(fl-^  !(■* 3  -0(0, 3  -  0$)  -  (1  -  2p)(5j<  +  hj  -m,!j  -  <t>jj  +  2c0,^3) 


-2v8ij(<t>,33  -  0,33  +  2c 0,333)  +  2oC30„y33] 


n(\  +  v)e 
2ir(l  -  v) 


U>,!j  +  4>.(j  -  2(1  -  2^(63,  +  53 j  -  1  )d>Jl  -  4i»8y0,33 


+  2jc30,"3). 


For  points  inside  fij,  the  elastic  stress  a*j  is  given  by 


(j*  =  a  -  a** 


(29) 


=  (o'j  -  <***)  +  a"  +  (T,-,  (30) 

where  -erf*  is  the  uniform  stress  that  exists  in  the  inclusion  caused  by  the  uniform  eigenstrain  ejj 
(Eq.  2).  The  stress  (o[j  -  of*)  is  the  uniform  stress  inside  the  inclusion  when  the  medium  is 
infinite.  The  solution  is  expressed  explicitly  by  Mura  ([13],  Eq.  11.20).  Equations  (5),  (12),  and 
(26)  give  the  stress  field  in  cylindrical  coordinates. 
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Seo  and  Mura’s  results  [12]  for  the  elastic  field  in  a  half  space  caused  by  an  ellipsoidal  inclusion 
with  uniform  dilatational  eigenstrain  (obtained  by  using  Mindlin’s  solution  [25]  for  Green’s  function 
in  the  half  space)  can  be  obtained  as  a  special  case  by  taking  b  =  0  (and  a]  =  a2)  in  Eqs.  (29)  and 
(30).  Mindlin  and  Cheng’s  results  [9]  for  a  sphere  can  also  be  obtained  as  a  special  case  by  taking 
a,  =  a2  =  a3  and  b  =  0  in  Eq.  (29). 


ELASTIC  STRAIN  ENERGY 

The  elastic  strain  energy  can  be  expressed  as 


W  =  —  —  (  o*  eT  dV 
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=  1 


where  u*  is  the  dilation  stress  field  in  the  inclusion.  It  is  given  bv 
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when  b  =  0,  the  strain  energy  obtained  is  the  same  as  that  obtained  by  Seo  and  Mura  [12]. 

THE  ELLIPSOIDAL  INHOMOGENEITY 

When  an  inhomogeneity  contains  an  eigenstrain,  it  is  called  an  inhomogeneity  inclusion. 
Eshelby  [6]  first  pointed  out  that  the  stress-field  changes  caused  by  an  inhomogeneity  when  the 
remotely  applied  stress  is  a,"  can  be  simulated  by  the  eigenstress  caused  by  an  inclusion,  if  the  eigen¬ 
strain  ejj  is  properly  chosen.  This  eigenstrain  is  sometimes  referred  to  as  the  equivalent  eigenstrain, 
or  the  equivalent  stress-free  transformation  strain.  For  a  given  uniformly  applied  stress  a/J  and  a  uni¬ 
form  eigenstrain  ej* ,  the  normal  components  of  the  equivalent  eigenstrains  e ?  are  given  by  [23] 


(X  -  X*)ec  +  \e‘ 


2 v-ejj  +  2(n*  -  n)  £ 
w  =  1 1 


Sijki  eh 


=  \*eT'  +  (X  -  \*)ea  +  2 n*e?  +  2 0*  -  (33) 


where  ij  =  11,  22,  33  and  kl  denotes  summation  over  11,  22,  33  only;  eT ,  e7*  and  ea  are  the  sum 
of  three  normal  components  of  strains  e ejj ,  and  e-j  respectively; 
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In  this  equation,  X  are  the  elastic  constants  of  the  matrix;  n* ,  X*  are  the  elastic  constants  of  the 
inhomogeneity;  and  /,,  /2,  lj,  and  Sljkt  are  constants  whose  values  depend  on  the  shape  of  the  inclu¬ 
sion  as  given  by  Eshelby  [6-8],  Some  detailed  expressions  for  these  constants  for  the  inclusions  of 
special  shapes  are  given  by  Mura  [13].  Therefore,  by  solving  the  set  of  three  simultaneous  equations 
in  Eq.  (33),  the  equivalent  eigenstrains  of ,  and  ^33  are  obtained  once  the  uniform  eigenstrain 
e/j'  and  uniformly  applied  stress  e°  are  given.  If  both  ej*  and  e-j  are  axisymmetric  for  an  axisym- 
metrical  inclusion,  die  resultant  equivalent  eigenstrain  ejj  is  also  axisymmetric  and  can  be  represented 
in  the  form  of  Eq.  (2).  Then  the  results  of  Eqs.  (3),  (5),  (12),  (26),  (27),  (29),  and  (30)  can  be 
applied  accordingly  to  solve  the  stress  field  and  strain  energy  of  an  axisymmetrical  inhomogeneous 
inclusion  in  the  half  space. 

SURFACE  DISTORTION  AND  DILATATION  FIELD 


The  roughness  of  solid  surfaces  is  a  second-order  effect,  but  it  has  profound  practical  conse¬ 
quences  in  many  fields  of  engineering  and  pure  science.  In  many  practical  situations,  the  presence  of 
inclusions  or  inhomogeneities  under  an  external  load  will  change  the  surface  profile.  The  displace¬ 
ment  of  the  free  surface  ( z  =  0)  solved  by  the  present  method  is: 


ur  = 
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The  presence  of  inclusions  or  inhomogeneities  under  an  external  load  will  also  produce  a  dilata- 
tional  field.  The  dilatational  field  in  the  matrix  obtained  in  the  present  study  is: 
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The  important  relationships  between  the  dilatation  field  and  the  equilibrium-concentration  distribution 
for  dilute  solutions  in  stressed  solid  arc  given  by  Li  [26]. 

SUMMARY 

The  stress  field  in  tne  half  space  (z  >  0)  caused  by  a  penny-shaped  inclusion  ft,  centered  at 
(0,  0,  c)  with  eigenstrain  =  6,,(e  +  bb (3)  is  found  by  the  superposition  of  the  following  three 
stress  fields:  (a)  the  stress  field  of  the  inclusion  ft,  centered  at  (0,0,c)  with  eigenstrain  ejj  in  an  infi¬ 
nite  medium;  (b)  the  stress  field  of  the  image  inclusion  ft2  centered  at  (0,  0,  -c)  with  eigenstrain 
—ejj  \  and  (c)  the  additional  fictitious  stress  field  that  makes  all  stress  fields  satisfy  the  equilibrium  and 
boundary  conditions. 

The  stress  field  of  the  penny-shaped  prismatic  inclusion  in  an  infinite  medium  obtained  by 
Eshelby  is  compared  with  the  stress  field  of  a  prismatic  loop  in  an  infinite  medium  as  obtained  by 
Kroupa  [21],  A  relationship  is  found  between  the  potential  function  </>  of  the  inclusion  and  the  integral 
function  /0_1 ,  which  involves  the  product  of  the  Bessel  functions  Jm  for  the  solution  of  the  prismatic 
loop. 
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The  fictitious  stress  field  is  solved  first  for  the  two-dimensional  problem  by  using  the  Hankel 
transformation  method  and  then  it  is  transformed  into  the  **  ree-dimensional  case  by  use  of  the  rela¬ 
tionship  between  <j>  and  /0~ 1 . 

The  solution  of  the  elastic  field  in  the  half  space  with  ellipsodial  inclusions  with  uniform  dilata- 
tional  eigenstrains  obtained  by  Seo  and  Mura  (1979)  has  been  rearranged  into  three  terms  correspond¬ 
ing  to  the  sress  field  of  the  inclusion  0,  in  an  infinite  medium  centered  at  (0,0, -c)  with  eigenstrain 
6tJ  e ,  the  stress  field  of  the  image  inclusion  fl2  centerd  at  (0,0, -c)  with  eigenstrain  -<5,,e ,  and  the 
additional  fictitious  stress  field.  It  has  also  been  shown  that  when  a{  =  a2,  Seo  and  Mura’s  results 
are  a  special  case  of  the  present  solution. 
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